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Consensus-based control for a network of 
diffusion PDEs with boundary local interaction 

Alessandro Pilloni, Alessandro Pisano, Yury Orlov and Elio Usai 


Abstract 

In this paper the problem of driving the state of a network of identical agents, modeled by boundary- 
controlled heat equations, towards a common steady-state profile is addressed. Decentralized consensus 
protocols are proposed to address two distinct problems. The first problem is that of steering the states 
of all agents towards the same constant steady-state profile which corresponds to the spatial average 
of the agents initial condition. A linear local interaction rule addressing this requirement is given. 
The second problem deals with the case where the controlled boundaries of the agents dynamics 
are corrupted by additive persistent disturbances. To achieve synchronization between agents, while 
completely rejecting the effect of the boundary disturbances, a nonlinear sliding-mode based consensus 
protocol is proposed. Performance of the proposed local interaction rules are analyzed by applying a 
Lyapunov-based approach. Simulation results are presented to support the effectiveness of the proposed 
algorithms. 


Index Terms 

Average consensus. Synchronization, Heat equation. Boundary control. Sliding-mode control. 


I. Introduction 

The problem of understanding when individual actions of interacting dynamical agents give 
rise to a coordinated collective behavior has received considerable attention in many research 
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fields. Examples can be found, e.g., in system biology [1], sensors networks [2], robotics [3], 
etc. 

In the systems and control literature, the usual mathematical setup underlying this class of 
problems refers to a group of agents, each one described by a dynamical system with one or more 
inputs, along with a communication network. Agents connected by a communication link are said 
to be neighbors, and can exchange information in either bidirectional or unidirectional manner. 
This raises the problem of designing decentralized local interaction policies (where each agent 
can only access neighbors information) in order to orchestrate the global coordinated behavior 
of the network. Within this framework, the consensus problem seeks to enforce agreement 
amongst the states of networked dynamical systems by penalizing their local disagreement with 
the neighboring nodes in a dynamic manner. The reader should refer, e.g., to [4], [5] for tutorial 
overviews of consensus-based control in the finite-dimensional setup. 

There are deep connections between the consensus problem and certain Partial Differential 
Equations (PDEs), the diffusion equation in particular [6]. Eor instance, discretizing in the spatial 
domain the one dimensional diffusion equation yields a high-dimensional system of networked 
first-order continuous-time integrators interacting through a linear Eaplacian-based consensus 
protocol. More generally, the application of finite-difference approximations of PDEs results 
in the spatial variables being mapped into the agent indexes, and the spatial derivatives being 
transformed into links between neighbors. Owing on the deep connections between the consensus 
problem and certain Partial Differential Equations (PDEs) (see e.g. [6]), some authors (see 
e-g- [7], [8]) have exploited (discretized forms of) several distributed parameter systems, such 
as advection and diffusion-advection equations, to derive more effective consensus protocols 
with improved convergence features. In spite of this intimate relationship between consensus 
algorithms and certain discretized PDEs, the consensus problem for a network of agents modeled 
as distributed parameter systems has not received yet the same level of attention than its finite¬ 
dimensional counterpart. 

The following papers [9], [10], [11], [12], [13], [14], [15], [16], which have investigated 
different aspects of consensus and synchronization in the distributed parameter systems setting, 
are worth to mention. In [9], [10], application of consensus to controlling mobile actuators in 
diffusion processes is discussed. The problem of designing consensus filters for state estimation 
in order to effectively integrate local information coming from a distributed spatial domain has 
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been widely addressed [11], [13], [14], both in the non-adaptive and adaptive setting. In [17], 
exact synchronization was achieved for a set of coupled wave processes, each one equipped with 
a boundary control input. 

In [18], [19], the authors present consensus algorithms in the framework of Multi-Agent 
Systems (MASs) modeled by continuous-time Partial difference Equations (PdEs) on graphs. 
Conceptually, this class of PdEs mimics PDEs in spatial domains having a graph structure, and 
many mathematical tools of functional analysis for studying PdEs are completely analogous to 
the ones developed for PDEs. In [20], [21], PdE-based model reference adaptive control laws 
are designed for a network of mobile agents to track desired deployment trajectories . 

Within a related framework, the recent work [15] studies the 3D agents deployment problem 
by PDE techniques, treating the agents as a continuum and modeling their interaction through a 
complex-valued diffusion-reaction PDE in a 2D spatial domain. An explicit backstepping-based 
local boundary control is designed to stabilize a variety of open-loop unstable deployment mani¬ 
folds. [22] presents a feedforward controller for multi-agent deployment by using a flatness-based 
motion planning method for PDEs. Reaction-advection-diffusion PDEs are used in [23] along 
with a backstopping design for leader-enabled agents deployment onto planar curves. Similarly, 
hyperbolic PDE models are used to design decentralized control laws for large vehicular platoons 
[24] and to analyze networks of oscillators [25]. 

In [16], synchronization and consensus problems have been studied for a network of agents 
modeled by a class of parabolic PDEs and communicating through undirected communication 
topologies. In [16], which appears to be more closely related to the present investigation among 
the existing references, some noticeable results are attained. Eirst, the case of all-to-all communi¬ 
cation between a set of identical agents is investigated and then the more realistic communication 
topology where each agents can only communicate with a limited set of neighboring agents 
is studied. In addition, a general abstract formulation of the underlying agents dynamics is 
introduced, and linear consensus controllers are presented to ensure agents agreement with a 
guaranteed convergence rate. It should be noted however that the common steady-state profile of 
the agents is not established, and furthermore no perturbations are allowed to affect the agents 
dynamics. 

The present work aims to address consensus and synchronization problems for a network of 
dynamical agents, communicating through an undirected and connected static topology, provided 
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that agents dynamics are governed by a class of diffusion PDEs with Neumann-type boundary 
actuation. The contribution of the paper is twofold. Firstly, a linear local interaction strategy 
whose implementation requires collocated boundary sensing only is proposed. With this strategy, 
it is shown that the agents states eventually converge “pointwise-in-space” towards a common 
constant distribution whose value is given by the spatial average of the agents’ initial conditions. 
Thus, the well-known average consensus algorithm is generalized from a network of integrators 
to the infinite-dimensional setting of networked heat processes. 

Secondly, the more complex scenario where the agents dynamics are perturbed by a class of 
boundary disturbances, is considered. Based on the second-order sliding-mode control approach 
[26], a nonlinear protocol is developed to extend the results of [27] from a network of double 
integrators to the infinite-dimensional framework of networked PDEs. A dynamic input extension, 
similar to that presented in [28] for stabilizing a unique perturbed diffusion PDF, results in 
continuous boundary control actions thereby alleviating chattering and yielding another step 
beyond [27]. It is demonstrated that the proposed nonlinear local interaction protocol, which only 
employs sensors located at the controlled boundaries, enforces the asymptotic synchronization 
between the agents states while rejecting the persistent matching boundary perturbations. 

The motivation to the present investigation comes, e.g., from networked systems of perturbed 
heat equations that can occur in modeling and controlling industrial furnaces. Heating of certain 
industrial furnaces (see, e.g., [29, Sect. l.A]) is made through electrically heated bars aiming 
to enforce a uniform temperature distribution inside the furnace. Considering these bars as a 
network of heaters and applying collaborative consensus-based synthesis might be useful in 
improving the overall performance of the furnaces. Exploiting the present results in specific 
application domains certainly requires additional work which is beyond the scope of the present 
paper. 

The paper is organized as follows. In Section II some mathematical preliminaries and useful 
properties and definitions are recalled. The linear average consensus algorithm is presented in 
Section III whereas the nonlinear algorithm, providing robust synchronization in the presence of 
boundary perturbations, is described in Section IV. Simulation results, supporting the proposed 
designs, are given in Section V, and conclusions and perspectives for next investigations are 
collected in the final Section VI. 
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II. Mathematical Preliminaries and Notations 
A. Useful definitions and properties 

The /p-norm and the /oo-norm of the real-valued N-dimensional veetor x= [ti, ... ,Tyv]^ G 
are defined as ||t||p = 1 < p < oo and ||a:||oo = maxi<KAf |.r 2 |, ■ ■ For 

the /i- and I 2 - norms, the following inequality holds [30]: 

||t||2 < ||t||i < sfN ||jc||2- (1) 


Let p, q > 1 be given, sueh that l/p-|-l/q = 1. Then the next ehain of inequalities is in foree 
[30] 


X y < ||t| 


lq< 


+ 


p q 

Operator sign(v), v G M, stands for the multi-valued funetion 


( 2 ) 


sign (v) G < 


1 

[- 1 , 1 ] 

-1 


if V > 0 

if V = 0 , 

if V < 0 


( 3 ) 


T 

whereas Sign (.r) stands for the veetor Sign (.r) = [sign {x\) , sign (X 2 ),..., sign (.ta^)] • 

The identity matrix of dimension N is denoted as InxN ^ whereas Ia? = [1,1, • • •, 1]^ G 

and Oa' = [0,0,...,0]^ G stand for the all-ones and all-zeros vectors. 

H''(0,1), with r = 0,1,2,..., denotes the Sobolev space of absolutely continuous scalar func¬ 
tions z{g) on the domain (0,1), with square integrable derivatives up to order i and the 

ff-norm ||z(-)llff(o,i) = \//o‘ELo f'*>(«)]" 

Then, the notations 


H''’^(0,1)= ^ 


H''(0,1) xH'"(0,l) X ...xH"(0,l) 


V 

N times 


( 4 ) 


and 


N 


|w(-)|| 






i=l 


2 

Hqo,!) 


( 5 ) 


for the corresponding norm of the vector w{g) = [wi{g),...,WN{g)Y ^ 1) are utilized. 

The simplified notation ||z(-)||h'' = lk(‘)Undo,!), IK(')IIh''A = IK(')IIh''A(o,i) will be adopted 
throughout. 
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For later use, an instrumental lemma is further presented: 

Lemma 1: Let b{g) E 1). Then, the following inequality holds: 

||fc(')ll?,o,»<2(||fc(0lli+l|fcs(-)l&,»). <' = 0,1 

Proof of Lemma 1: It was proven [28, Lemma 1] that, with reference to a generic scalar function 
z{g) E 1), the next estimate holds: 

lk(-)ll?io<2(z(/f + ||z,(.)||2o), / = 0,1. (6) 

Now let b{q) = [bi{g)M{‘;):---:bN{q)Y and b^{g) = [b^i{g),b^2{g),...,b^N{<;)Y where 
bk{g) £ H^(0,1) 'ik= By applying definition (5), the following chain of relations 

is derived by virtue of (6) specified with z(-) =bk{-): 

N N 

IIM0llHO,- = III^K0llHO<2L(^K0" + II^^K0llHo)=2(IIM0lli + ll^a0llHO-)^ ^- = 0,1. 

f=i j=i 

(7) 

Lemma 1 is proved. □ 

B. Algebraic Graph Theory definitions and properties 

We consider a set of N dynamical agents along with an undirected static communication topol¬ 
ogy represented by the graph Q{V,S), where V = {1,... ,N} is the set of vertices representing 
agents and £ C {V x V} is the set of edges representing the information flow among the agents. 
J\fi = {j eV (ij) e£} denotes the set of neighbors of agent i. The topological structure of ^ 
is encoded in the so-called Laplacian Matrix C = fij] E where 

\M\ if i = i 

iij := < —1 if (i,j)ES (8) 

0 otherwise 

For undirected connected graphs, the matrix C is symmetric and positive semi-definite [5], the 
properties 

CIn = C^In = ^n ( 9 ) 

hold by construction, and the corresponding eigenvalues A,-, i E V, are such that 0 = Ai < A 2 < 
...< The smallest nonzero eigenvalue A 2 is known as algebraic connectivity of Next 
lemma presents useful properties of vector norms involving the Laplacian matrix of the graph. 
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Lemma 2: For an undirected connected graph with Laplacian matrix C, and with reference 
to any vector x G such that \^x = 0, the next relations are in force 




AtvI |;c| I2 > x^Cx > A2I | 7 c| I2 

,,2 ^ ^^2 ll-^lli 

Cx\\i > XiWAh 


\l>\\Cx\\^,>XA 


( 10 ) 

( 11 ) 

( 12 ) 


Proof of Lemma 2: The left inequality in (10) comes from well-known properties of quadratic 
norms. The right inequality in (10) was proven in [5, Th. 3]. To derive (11), observe that ||jC.x :||2 = 
\/x^C^x. The eigenvalues {O, Af, ..., of are straightforwardly derived by squaring 
those of C. Thus, the left inequality of (11) follows from well-known properties of quadratic 
norms. Additionally, is symmetric and such that CAXjx = Oa^, thus Ia^ is the eigenvector 
associated to the zero eigenvalue of CA. Therefore, the right inequality of (11) follows from the 
Courant-Fisher Theorem that can be found, e.g., in [31]. To reproduce (12), it suffices to conclude 
from (1) that ||jCa3||i > ||;Ca7||2 and then, by applying (11), to derive that ||jCa3||2 > A 2 ||a?|| 2 - 
Lemma 2 is proved. □ 


III. Average consensus for networked heat processes 


A network of N dynamical agents whose communication topology is described by an undi¬ 
rected connected static graph ^(V, £) is under study. The f-th agent has state Qfgp), i G V, with 
the spatial variable G (0,1) and time variable t > 0. Let e(g,o = ei(g,o, e2(?T), • ■e^(?T) 
be the vector collecting the states of all agents, and let the dynamics of Q{(^d) be governed by 
the vector heat equation 

Qr{gd) = e-Q^Aif)^ ( 13 ) 


The scalar parameter 0 GM+ is a positive unknown coefficient, called “diffusivity parameter”, 
which is supposed to be identical for all agents. Throughout, Neumann-type Boundary Conditions 
(BCs) of the form 


are considered, where U(t) 
ary vector control input). 


2 ^ 0,0 = 0 , QAht) = Uit), 


(14) 


T 


Ui{t),U2{t),...,UN{t) 


G is a modifiable source term (bound- 
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The Initial Conditions (ICs) are 


Qig,0) = Qoig) (15) 

To deal with classical solutions of class 1), the admissible initial functions are specified 

by the next assumption. 

Assumption 1: The initial function Qo{g) in the ICs (15) is assumed to be of class H^’^(0,1) 
and compatible to the BCs 2og(0) = 0 2og(l) = U{f)). 

The objective of the present section is to introduce a linear local interaction strategy providing 
closed-loop stability and the point-wise consensus condition 

lime(g,0 = e*-lyv, Vge(0,l), (16) 

r—>oo 

where the constant 

e* = -f /' i«eo(s)ds (17) 

N Jo 

corresponds to the spatial averaging of the agents initial conditions. 

To achieve the control goal, the local interaction protocol 

U{t) = -CQ{l,t) (18) 

is proposed. Under the assumptions, imposed on the ICs and BCs, the well-posedness of the 
system in question is straightforwardly verified by applying [32, Theorem 2.1.10] to the classical 
solutions of the homogeneous linear Boundary-Value Problem (BVP) (13)-(15), (18). 

We are now in a position to state the first main result of this paper. 

Theorem 1: Consider the multi-agent system (13)-(15), with Assumption 1, communicating 
through an undirected connected static graph with Laplacian matrix C. Let it be subject to the 
boundary local interaction control strategy (18). Then, the closed-loop system is stable in the 
space //^(0,1) and the average consensus condition (16)-(17) is achieved. □ 

Proof of Theorem 1: The stability of the closed-loop BVP is established by involving the 
Lyapunov function 

Viit)=l- Q\^d)Q{^d)d^ ( 19 ) 

2 Jo 
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whose time derivative, estimated along the solutions of the BVP (13)-(15), (18), is non-positive 
definite: 

Mt) = e f = -e\\Q^{-,t)\\lm-QQ{uYCQ{\,t) < 

J 0 

( 20 ) 

To get (20) integration by parts, and BCs (14),(18), and the semi-definite positiveness of the 
Laplacian matrix C were utilized. Relation (20) ensures that the system is stable in the space 
//®(0,1). Since the BVP is linear, and the ICs (15) are of class //^(0,1), the stability remains 
in force in the space //^(0,1). 

Next, let us note that the eigenspace of the closed-loop BVP (13)-(15), (18), associated with 
the zero eigenvalue, is one dimensional and it is spanned by the uniform distribution Q{g) = 1a?. 
In turn, all the remaining eigenvalues are strictly negative because the BVP (13)-(15), (18) has 
been shown to be stable. 

Furthermore, one observes that the projection 

2 ( 21 ) 

N Jo 

of the solution of the closed-loop system to the eigenspace, associated to the zero eigenvalue, 
remains constant, whereas all the remaining modes tend to zero because the other eigenvalues are 
strictly negative. It follows that the state eventually converges point-wise to the constant 

spatial distribution 

i f\l,Qoi^)d^lN = Q*-lN (22) 

N Jo 

thereby establishing relations (16)-(17). This completes the proof of Theorem 1. □. 


IV. Robust synchronization for networked heat processes with perturbations 

A perturbed version of the BVP (I3)-(I5), with the only difference in the BCs (14) which 
now take the perturbed form 


25(0,0 = 0, 25(1,0 = t/(0+'i^(0, 


(23) 


is under investigation, where ^^(0 = 


pA? 


represents an uncertain. 


sufficiently smooth, persistent disturbance. 

The class of admissible ICs and disturbances is specified by the next assumption. 
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Assumption 2: The initial function Qo{(^) is assumed to be of class 1) and compatible 

to the perturbed BCs 20^(0) = 0? QoqiX) = whereas the disturbance is supposed to 

be twice continuously differentiable, and there exists an a-priori known constant IT > 0 such that 

n. (24) 

Note that for technical reasons (see Remark 1 below) higher degree of smoothness of the ICs 
is required in the present perturbed scenario. The objective of the present section is to develop 
a local interaction strategy providing the attainment of the synchronization condition 

lim \Qfqp)-Qj{<;p)\ = 0, V i,j e V,V$ e (0,1), (25) 

despite the presence of the uncertain boundary disturbance 'T(t) of arbitrary shape and possibly 
unbounded in magnitude. 

To achieve the control goal, the following dynamic local interaction protocol 


U{t) = Ui{t) + U2{t) (26) 

is proposed, with 

f/i (t) = - aSign (£e( 1 , t)) - hSign (£a ( 1, t)) (27) 

U2it)=-WrCQ{ft)-W2-CQ,ift)-W3-Qtift) (28) 

t/l(0)=t/2(0)=0yv. (29) 


The initial values f/i( 0 ),t/ 2 ( 0 ) are all set to zero to verify the compatibility^ Qogi^) = U{0) + 
'T(O) to the BCs (23) at g = 1. In (27)-(28), a, b, Wi, W 2 and IT 3 are nonnegative tuning constants 
subject to certain design inequalities that will be constructively derived in the sequel. 

It is worth to note that the discontinuities affect the time derivative of the boundary control 
vector, whereas the boundary control signal is smoothed by passing these discontinuities through 
an integrator, thereby alleviating chattering. 

Remark 1: Although the state derivative is normally not permitted to be employed in the 
synthesis (as it generally induces algebraic loops), its use becomes acceptable when dynamic 
input extension is performed, what is indeed the case of the present dynamic synthesis where 
the input signal passes through an integrator. By virtue of this, the system state is augmented by 

^See, e.g., [33] for the need of certain compatibility conditions in the dynamic boundary control synthesis. 
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Qt being viewed as a component of the augmented state vector (2, Qt) G 1) x 1) 

which is particularly why the initial function Qo{(^) was assumed to be of class 1). □ 

The well-posedness of the underlying closed-loop system, under the assumptions, imposed 
on the ICs and BCs, is actually verifiable in accordance with [32, Theorem 3.3.3] by taking 
into account that the dynamic local interaction rule (26)-(28) is twice piece-wise continuously 
differentiable along the state trajectories. Thus, in the remainder, it is assumed the following: 

Assumption 3: The closed loop networked system (13)-(15), (23), (26)-(29) possesses a 
unique Filippov solution 2(-,t) G 1) and its time derivative = Qt{-,t) ^ 1) 

verifies the auxiliary boundary-value problem 



(30) 

2^0,0 = 0, Z^{l,t) = U{t)+^it), 

(31) 

Z(g,O) = 02og5(g)GH4^(O,l). 

(32) 


Extension of the Filippov solution concept towards the infinite dimensional setting can be 
found, e.g., in [34]. Notice that (30)-(31) are formally obtained by differentiating (13)-(15), 
(23), in the time variable t, whereas the IC (32) is straightforwardly derived from (13) and (15). 

It is customary [5] to formalize the achievement of consensus through the annihilation of 
appropriate (N-dimensional) “disagreement” vectors. Generalizing [16], the following distributed 


disagreement vectors 

=-Cc Qi'd), (33) 

d),■■-1 {■ d)V = Shi-d) = die Qti'd), (34) 

die = (^NxN -, (35) 

will be considered in the present investigation for analysis purposes. The next properties hold 
due to (9) and (33)-(35) 

llSiigp) = ll52igd) =0, Vg G [0,1], (36) 

CCe = dledl = dl (37) 
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thereby implying that 


Cdi{g,t) = CCeQ{g,t) = CQ{g,t), 
C&2{g,t) = CCcQt{g,t) = CQt{g,t) 


The BVP governing the dynamics of the disagreement vectors now reads as 

du{g,t) = &i{g,t), 

&it{g,t) = ed2^g{g,t), 

52c(0,t) =0 

^ (41) 

&2g{ht) = Ce 

5i(g,o) = /:ceo($)GH4-^(o,i) 

&i{g,0) = eCcQo^^{g)en^^N{0,\) 

Presenting the second main result of this paper is preceded by the following instrumental 
lemma. 

Lemma 3: The functional 

V'(5,,52) = e<.||£i)i(l,<)lli + leM'ilK5|(l.<)ll2 + ^ f (43) 

2 2 ./o 

being computed on the solutions (5i(-,t), 52(-,t)) of the BVP (40)-(42), is equivalent to the 
h 2’^(0, 1) X 1) norm of these solutions in the sense that 

/ N 


for an arbitrary solution (5i (•, t), of (40)-(42), for all t > 0, and for some positive constants 

7)1 and r\ 2 . 

Proof of Lemma 3: It is preliminarily demonstrated that the condition 

ai -^(51,52) < y(5i, 52) < a2 -1^(51,52) (45) 


holds, where ai and a 2 are positive constants and 

V'((),,52) = e<.||5,(l,<)l|i + leM'il|5i(l.<)ll2 + ^l|52(-.<)llH.,»- 
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By considering (12) and the second inequality of (1), both specialized with 3c = 5i(l,t), one 
derives 

^.||5i(l.<)|l, < l|C5,(l.<)lli < l|£||il|5i(l.<)lli (46) 

Specializing (11) with ;c = 5i(l,t), and (10) with 3c = one obtains 

Aj2||5,(l,,)||Bll4:5i(l,<)lli<K«ll«i(l.<)ll2 (47) 

(48) 

Noticing that, by construction, Jq \\ 52 {^,t)\\\d^ = ||52(-,t)||^o.)V2 the next estimate is derived 
after spatial integration of all terms in (48) 

32||52(-,»)IIh».»< <3A.||52(-.»)llft«.»- (49) 

By (46)-(47) and (49), relation (45) is derived with the positive constants Cti = min{A 2 / \/N, A|} 
and a 2 = max{||/^||i , Aa?}. Furthermore, by (1) and (5), functional y(5i, 62 ) can be rewritten 

as follows: 

N 

y(5i,52)=£v^-(5i,-,52,-) (50) 

i=l 

where 

y = 0a|5i,-(l,OI + ^0W^i5n(l,t)" + ^||52/(-,t)||2o. 

From that, by applying [28, Lemma 2], the next estimate 

)3|(l|5ii(^0ll?,2+l|52i(^0llH»)<>i<ft(l|5ii('.<)llH2+IIS2/(',<)ll?,o+ll«i/(',<)llH2) (51) 

holds for some positive constants j8i and ^ 2 - Then, by combining (50) and (51), and by applying 
definition (5), it results 

)3l(l|5l(^0ll^2,» + l|52(^»)||ft0,,)<V(«l.S2)<)52M|«l(^»)||^,2,»+||S2(^»)||^,,,, + |;il«l.('.»)llH4 

(52) 

Finally, relation (44), being specified with the constants rji = aiPi and r ]2 = CC 2 P 2 , is straight¬ 
forwardly derived by combining (45), (50) and (52). Lemma 3 is proved. □ 


Next theorem presents the second main result of this work. 
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Theorem 2: Consider the perturbed multi-agent system (13), (15), (23), with Assumptions 
2 and 3, communicating through an undirected connected static graph with Laplacian matrix 
C. Let the boundary local interaction strategy (26)-(28) be applied, with the tuning parameters 
selected according to 

a>b + Yl, b>Yl, ILi > 0, 1^2 > 0, IL? > 0 (53) 

Then, condition (25) is achieved. □ 

Proof of Theorem 2: See Appendix for the proof. 

Remark 2: In contrast to the average consensus result, outlined in Section II, in this case the 
steady-state common profile reached by the agents cannot be predicted a-priori and it turns out 
to depend not only on the agents initial conditions but also on the actual controller parameters 
and disturbance vector. For this reason, the term “robust synchronization” has been adopted in 
the present scenario to describe the underlying result, as opposed to the word ’’consensus” that 
mostly refers, in the literature, to situations where the steady state behaviour of the agents is 
determined a priori and the problem is that of enforcing it in a decentralized manner. □ 

V. Simulation Results 

In the present section, simulation results are presented to illustrate the performance of the 
proposed protocols. The connected network of A = 10 agents displayed in Figure 1 is considered, 
with the diffusivity parameter 0 = 1. 



Fig. 1. The considered network topology. 
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Fig. 2. TEST 1; Spatiotemporal profiles (left) and 2 io(iJT) (right). 

For solving the resulting closed-loop system of coupled PDFs, the spatial domain g G [0,1] 
has been discretized by the standard finite-difference approximation method considering n = 30 
uniformly spaced solution nodes. The resulting finite dimensional system of coupled ODEs is 
then solved by means of the Euler fixed-step solver with sampling-step Ts — 10^"^. 

A. Average consensus 

System (13)-(14), coupled with the local interaction protocol (18), is under investigation. In the 
first simulation run (TEST 1), spatially varying ICs have been selected as follows: Qi{g,0) = 

10 + (Oicos(3;rg), 22 ($, 0 ) = 10 + © 2 Cos( 3 ;rg), ^ 3 (^, 0 ) = 8 + (O 3 cos( 3 ;rg), 24($,0) = 10 + 
m 4 Cos( 3 ;r<g), 25(<g,0) =6 + 0 ) 5 cos(3;rg), 26(<;,0) = 10+ci)6cos(3;rg), Qi{g,Q) = 10 +ci) 7 cos( 3 ;rg), 
28($,0) = -5 + ( 08 Cos( 3 ;rg), Q^{g,Qi) = 10 + (O 9 cos( 3 ;rg), 2io(g,0) = 10 + (Oiocos(2.5;rg), 
with G)/ = 1 +4(/ — l)/9, i= 1,2..., 10. The corresponding spatial average Q*, evaluated ac¬ 
cording to (17), is Q* = 1.9631. Figure 2 shows the spatiotemporal evolutions of the states 
Qe{<i,t) and 2 io($T)- The steady-state profile of both agents is constant and takes the expected 
pre-computed value Q*. Figure 3 displays the time evolution of the disagreement vector norm 

11 5i (•, t) ||h 2 ,io, which tends to zero thereby confirming that all agents reach a common steady-state 
profile. 

B. Robust Synchronization 

The performance of the robust synchronization protocol (26)-(28), presented in Section IV, 
is now verified with reference to the perturbed PDFs (13), (15), (23). The entries ^\fi{t) of the 


September 24, 2015 


DRAFT 




































16 



Time [sec] 


Fig. 3. TEST 1: Time evolution of the disagreement vector norm ||5i (•,t)||jj 2 .'v- 

disturbance vector 'P(t) are selected as !//,■(?) = Aki-t-\-sm{kiKt) (i = 1 , 2 , where the 

coefficients ki are randomly chosen in the interval [0,2]. The considered disturbance, which is 
unbounded in magnitude as time grows, meets the restriction (24) with a constant upper-bound 
constant n = 2;r + 8 . The chosen ICs are 

Qi{g, 0) = 10 + {i - 4.5) cos(4;r$) / = 1,2,..., 10. 

The tuning parameters a = 40,b = 20, ITi = W 2 = IT 3 = 5, were chosen according to Theorem 2. 
This simulation run is referred to as TEST 2. Figure 4 shows the spatiotemporal evolution of 
the state 26(<JT) and of the state mismatch Qioig^t) — It is clear that both states 

converge towards the same steady-state profile. The time evolution of the disagreement vector 
norm || 5 i(-,t)||j 32 ,io, shown in Figure 5-left, tends to zero as shown in the Theorem 2. To verify 
the conservativeness of the approach, another simulation run, called TEST 3, was made where the 
terms i = 1 , 2 ,..., 10 (a,- being randomly chosen constants in the interval [ 0 , 20 ]) has been 
added to the disturbance entries V/,(t) used in TEST 2. Due to the insertion of these additional 
terms, not only the magnitudes of the disturbances are increasing with time but also those of 
their derivatives are. Thus, the tuning conditions (53) are deliberately violated. Figure5-right 
depicts the resulting divergence trend of the disagreement vector norm ||5i(-,t)||H2.io, thereby 
supporting the theoretical results. 


VI. Conclusion 

In this work, an infinite-dimensional counterpart of the well-known finite-dimensional average 
consensus algorithm has been derived with reference to an unperturbed network of diffusion 
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Fig. 4. TEST 2: Spatiotemporal profiles (left) and 26(ST) ~ 2lo(irT) (right). 




Fig. 5. Time evolution of the norm ||5i(-,f)|lH2.'' ir* TEST 2 (left) and TEST 3 (right). 

processes under a linear deeentralized loeal interaetion poliey. Along with this, the problem 
of guaranteeing the asymptotie agreement among the agents’ states while rejeeting a elass of 
persistent and possibly unbounded disturbanees has been taekled by devising a nonlinear loeal 
interaetion poliey based on the seeond-order sliding-mode eontrol approaeh. 

Future aetivities will be targeted to relaxing the topologieal restrietions on the network strueture 
by eovering, e.g., direeted and possibly switehing eommunieation graphs. Further investigation 
is ealled for the extension of these results to more general elasses of (possibly non identieal) dis¬ 
tributed parameters agents’ dynamies and additionally to more general eonsensus-based problems 
in the infinite dimensional setting sueh as e.g. leader following. 
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Appendix 

Proof of Theorem 2 

Consider the Lyapunov function (43). Its time derivative is 

V{t)=ea& 2 {l:tfCSign{C5,{\,t)) + eWMl:tf[\i^jfCd 2 t{^j)d^ (54) 

Jo 

Substituting (40) into the last term of (54), and performing integration by parts in light of 
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(41), yield 

= e5,{i,tfc52i{i,t) -ej^ S2t(i,tf‘^hi{i,t)di (55) 

Owing on the spatial differentiation of equation (36), relation (10) specified with x= 82 ^ holds 
true. Thus, the last integral term of (55) can be estimated as follows 

-e£d2^(S,tVcd2f(^,i)d^ <-eX2\\S2i{;i)fa,,« (56) 

Substituting the BCs (41) and the controller equations (26)-(28) into the right-hand side of 
(55) one derives 

052 ( 1 , 0 ^^ 525 ( 1,0 =- 0 a- 62 ( 1 ,O^^Sign(£ 5 i(l,t))- 0 Z 7 ||£ 52 (l,Olli- 0 W^i 52 (l,O^^^ 5 i(l,t) 

-01T2||£52(l,Oll2-01^352(l,O^^a(l,O+ 052(1,O^^W (57) 

Employing (39), evaluated at g = 1, and considering (10) specified with jc = 52(1,0, one gets 

-01T352(l,O^^a(l,O < - 0 I+ 3 A 2 ■ ||52(l,0lli (58) 

By applying (2), specified with .r = '^(O, y = /^52(1,0, P = °°, ^ = 1, and taking into account 
(24), the magnitude of the last term in the right hand side of (57) is estimated as follows 

1052(1,0^^^1 = 01'i^^(0^52(1,01 < 0n||)C52(l,Olli (59) 

Substituting (55)-(59) into (54), and considering relation A|||5i(l,0ll2 + Il^5i(l,0ll2, which 
is verified by applying (11) specialized with .r = 52(1,0, the next inequality 

V(O<-0-(5-n)-||)C52(l,Olli-0A2-||525(-,OllHo.iv 

-01T2A|-1152(1,0112-01+3^2-1152(1,0112 (60) 

is concluded after some straightforward manipulations. Due to (60) and (53), the time derivative 
of the Lyapunov functional V{t), being computed along the solutions of the closed-loop system, 
is negative semi-definite, and V{t) is thus a non-increasing function of time. Thereby, the set 

= {(5i, 52) G X : y (5 i, 52) < R} (61) 
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specified for an arbitrary ^ > y(0), is invariant. By exploiting the invariance of the domain 
the next estimates are derived by straightforward manipulations of the inequality y(-) <R in 
light of (43), (47) and (49): 


l|£'Si(l.<)lli < 

R/6a 

(62) 

ll'Si(i.<)lli < 

2R/eWiXi 

(63) 

1 <52(-,t)| ^o,iV < 

2RIX2 

(64) 

Now consider the “augmented” functional 



VRit)=Vit) + KR-V{t) 


(65) 

v{t) = ^ew2-\\c5i{i,t)\\l+ 

[' d,il,tVCd2{^d)d^ 

Jo 

(66) 


where Kr is a sufficiently small positive constant to subsequently be specified. The next estimation 


holds 

^0W2-|K5i(l,t)||2>l0iy2A|-||5i(l,O|li (67) 


whereas by (1), (2) and (62), the second term in the right-hand side of (66) is manipulated as 

Considering (67) and (68) along with (65), (43), yield 

£5|(l.<)lll + ^^^l|S2(',<)llH0.» + 5e)(|(W'i + li-KW'2)||5i(l,<)ll2- (69) 

Thus, the positive definitiveness of V/?(t) is guaranteed by selecting the positive constant Kr 
small enough according to 

Kr < min|2^^^2} 


VR{t) > 


9a 


KrR 

26a 


. In particular, in the invariant domain the augmented functional VR{t) turns out to be lower 
estimated in terms of V (t): 

Oa A 2 Wi 




0a-^’^2-KR’Wi + KRW2 


V(l). 


(70) 


Differentiating (65) along the solutions of (40)-(42), and exploiting the identity jCCc = jC in 
(37), it yields 

VR(t) = V{t) + KRBW2&2{i:t)C^d,{\,t) + Kr [' 0 5i (1, t) ^£^2,^^ ((^ , + Kr [' 52ihtfC52i^j)d^. 


(71) 
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Employing (26)-(28) and the BCs (41) one gets 

£05i(l,t)^/:52,^^(<^,t)J<^ = 05i(l,t)^/:52,5(1,0 =-0a-||Z:5i(l,OI|i-0(^-5i(l,O/:Sign(/:52(l,O) 
- 01Ei-||£5i(l,Olli-01E2-5i(l,O^/:" 52(1,0-0^5-51(1,0^^52(1,0+ 0-51(1, 

(72) 

Then, utilizing (2) and (63), the sign-indefinite terms in (72) are estimated as follows 


|-5i(l,0/:Sign(/:52(l,0)|<IK5i(l,0l|i , 
-5i(l,0^r52(l,0l<l|5i(l,0ll2lK52(l,0ll2 


<||5i(1,0I|2|K52(1,0IIi< 


2R 


01 TiA| 


■IK52(i,0IIi , 


5i(i,o^/:WI<ll'i^(OII4^5i(i,i 


(73) 


(74) 


(75) 


By (64) and the Holder integral inequality, the last integral term in the right hand side of (71) 
is estimated as 


70 




<,/^-||£52(1.<)IIi 


d2(\,ifcd2(i,t)di 

By (71)-(76), and considering (11) and (60), one manipulates (71) as follows: 

Vff(0 < -0 (b-n-KR^lR/e^X^ ||£52(l,/)lli + K-ffW30^2/?/0WiA|.||£52(l,Olli 


(76) 


-K-«0 (a-fe-n). ||£5l(l,0||l-0^2-||52,s(-,f)llHOA-0(lE2A| + H6A2 )-1152(1,0112-1151(1,0112 

(77) 


By Lemma 1, specialized for 5(-) = 52 (-) and / = 1, the next estimate 

-0(lT2A| + lT3A2)||52(l,Olli-0A2||525(-,Oll2<-C4||52(-,OllHO,^ (78) 

is obtained with 

C 4 = 0 A 2 min { 1 , (IT 2 A 2 + IT 3 )} 

. By substituting (78) into the right hand side of (77), the inequality 

VR{t) <-ci • ||/:5 i(1,0I|i-C2-|K52(1,0IIi-C 3- ||5i(l,0ll2-C4- ||52(-,0 IIho.a' (79) 

is finally obtained, with the coefficients 

Cl = KrO {a — b— Yi), C3 = k-^OITiAI, 
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—(-"-(vS-se)) 

It is clear that due to the proposed specifications of constants ci,C 2 ,C 3 , all terms, appearing in 
the right-hand side of (79), are nonpositive provided that the tuning conditions (53), imposed 
on the controller parameters, hold and the next more restrictive condition 

T b-n 

^ I 2R I ^3 

^2 V 


is additionally satisfied. It then follows from (79) that 


vi(<)<-ri(llc«i(i.<)lli + ll«i(i.<)lli+ll«2(-,<)ll?,o,») (80) 


with 7 i = min{ci,C 3 ,C 4 } > 0. By (43) combined with the first inequalities of (10) and (11) it 
yields 

y(O<0a||r5i(l,t)||i-f^0WiA^||5i(l,Olli + ^||52(-,t)||HO.(v, 

whereas by (10), (62) and (1), along with property (2) specialized with v = Cdi(lj) and y = 
52(iJ,t), one derives that 

Finally, substituting the last two estimations in (65) one obtains 


Vr (0 < 72(||£5i(l,t)||i -t- ||5i(l,t)||2+ ||52(-,t)||po(Q i)]«) 


(81) 


where 


■ (a + ?Q;VF2)., ^ ^ 

72 = mm{0a - ,---, ---} > 0. 


20a’ 2 

Thus, one derives from (80) and (81) that 


VR{t) < -pR-VR{t), PR = YiIY 2, 

thereby concluding the exponential decay of VR{t), initialized within the invariant set in 
(61). 

To complete the proof, it remains to note that due to the upper estimate (70) the functional V(t) 
decays, too. By applying Lemma 3, the local asymptotic stability of (40)-(42) is then established 
in the space //^’^(0,1) x//®’^(0,1) for the initial set (61). Since (61) can be specified with 
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an arbitrarily large R > 0, thus capturing an arbitrarily large initial domain, and the tuning 
conditions (53) do not depend on R, the global asymptotic stability is then concluded in the 
space 1) xH®’^(0,1). It follows from (44) that ||5i(-,t)||H2,Ar asymptotically vanishes 

too, which results in the following component-wise relations 

lim||5i,-(-,0||H2 = 0, V/eV, (82) 

r—>oo 

It is well known [35] that the Sobolev space //^(0,1) is continuously embedded in the Banach 
space C(0,1) equipped with the supremum norm. In other words, there exists a constant M > 0 
such that 

sup^e[0,i]\Sii{^,t)\<M\\5u{-,t)\\^2, V/gV, (83) 

Thus, one concludes the spatially point-wise decay of all entries of (5i (-,?). This property, cou¬ 
pled to the identity Qi{g,t) — Qj{g,t) = 5ii{g,t) — 5ij{g,t), yields (25). The proof of Theorem 2 
is completed. □ 
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